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HIGHLIGHTS 


•  We  present  two  energy  supervisory  controls  for  improving  long-term  operation  of  off-grid  hybrid  systems. 

•  Battery  SOC,  hydrogen  tank  level  and  remaining  energy  of  storage  devices  were  considered  by  the  controls. 

•  The  controls  were  tested  by  simulation  along  the  expected  lifetime  of  the  hybrid  system. 

•  The  results  demonstrated  the  feasibility  of  the  controls  for  off-grid  applications. 
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This  paper  presents  two  novel  hourly  energy  supervisory  controls  (ESC)  for  improving  long-term  oper¬ 
ation  of  off-grid  hybrid  systems  (HS)  integrating  renewable  energy  sources  (wind  turbine  and  photo¬ 
voltaic  solar  panels),  hydrogen  system  (fuel  cell,  hydrogen  tank  and  electrolyzer)  and  battery.  The  first 
ESC  tries  to  improve  the  power  supplied  by  the  HS  and  the  power  stored  in  the  battery  and/or  in  the 
hydrogen  tank,  whereas  the  second  one  tries  to  minimize  the  number  of  needed  elements  (batteries,  fuel 
cells  and  electrolyzers)  throughout  the  expected  life  of  the  HS  (25  years).  Moreover,  in  both  ESC,  the 
battery  state-of-charge  (SOC)  and  the  hydrogen  tank  level  are  controlled  and  maintained  between  op¬ 
timum  operating  margins.  Finally,  a  comparative  study  between  the  controls  is  carried  out  by  models  of 
the  commercially  available  components  used  in  the  HS  under  study  in  this  work.  These  ESC  are  also 
compared  with  a  third  ESC,  already  published  by  the  authors,  and  based  on  reducing  the  utilization  costs 
of  the  energy  storage  devices.  The  comparative  study  proves  the  right  performance  of  the  ESC  and  their 
differences. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nowadays  problems  related  with  the  lack  of  electrical  grid  in 
inaccessible  areas  due  to  the  operating  cost  can  be  solved  installing 
off-grid  system.  Usually,  since  years  ago,  most  of  off-grid  systems 
have  been  based  on  diesel  engines  as  primary  energy  sources  [1,2]. 
However,  due  to  the  low  maintenance  and  operation  cost  and  the 
long  expected  life,  off-grid  system  based  on  renewable  energy  is 
becoming  an  alternative  to  the  classical  off-grid  system.  Normally, 
in  these  new  systems,  the  renewable  energy  is  composed  of 


*  Corresponding  author.  Tel.:  +34  956  02  81  66;  fax:  +34  956  02  80  01. 

E-mail  addresses:  pablo.garcia@uca.es  (R  Garcia),  jtorregl@ujaen.es  (J. 
R  Torreglosa),  luis.fernandez@uca.es  (L.M.  Fernandez),  fjurado@ujaen.es  (F.  Jurado). 


photovoltaic  solar  panels  (PV),  wind  turbines  (WT)  or  both,  so  that 
they  do  not  emit  greenhouse  gases,  they  are  low  pollution  energy 
sources  [3]  and  they  avoid  the  use  of  fossil  fuel  4]. 

However,  the  energy  production  from  the  sun  and  from  the 
wind  do  not  always  coincide  with  the  power  demand,  and  there¬ 
fore,  energy  secondary  systems  must  be  used  together  with  the 
renewable  energy  sources  for  energy  support/storage  [5,6  .  These 
secondary  systems  are  typically  batteries  and  hydrogen  systems, 
composed  by  fuel  cell  (FC)  electrolyzer  and  hydrogen  tank,  but  also 
super-capacitors  or  superconducting  magnetic  energy  storage 
(SMES)  [7]  can  be  used  [8  .  For  example,  the  storage  and  support 
system  was  integrated  by  hydrogen  system  (FC,  electrolyzer  and 
hydrogen  tank)  in  Ref.  [9];  by  battery  pack  in  Ref.  [10  ;  and  by 
hydrogen  system  and  battery  in  Ref.  [11].  Otherwise,  some  studies 
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Nomenclature 

Niz 

number  of  electrolyzers  used  throughout  the  lifetime 

CAPH2 

tank  nominal  capacity,  (kg) 

P* 

1  bat 

of  the  HS,  (-) 

ideal  battery  power,  (W) 

DODnom 

nominal  depth  of  discharge  for  the  battery,  (%) 

p* 

1  fc 

ideal  fuel  cell  power,  (W) 

^bat.char 

battery  charge  energy,  (Wh) 

p* 

1  lz 

ideal  electrolyzer  power,  (W) 

^bat.dis 

battery  discharge  energy,  (Wh) 

Pi 

time  that  the  ECS  have  been  working  in  the  state  1,  (%) 

cnom 

cbat 

battery  nominal  capacity,  (Wh) 

P48 

time  that  the  ECS  have  been  working  in  the  states  4  and 

Ef 

c 

FC  energy,  (Wh) 

8,  (%) 

E\ 

ow,H2 

hydrogen  lower  heating  value,  (Wh  kg-1) 

Ps 

time  that  the  ECS  have  been  working  in  the  state  5,  (%) 

E\ 

z 

electrolyzer  energy,  (Wh) 

pav 

bat,charg 

battery  available  power  during  a  charge  in  the  next 

ESC 

energy  supervisory  control 

hour,  (W) 

ESC,1 

superscript  for  the  equations  related  to  the  first  ESC 

pav 

bat,dis 

battery  available  power  during  a  discharge  in  the  next 

ESC, 2 

superscript  for  the  equations  related  to  the  second  ESC 

hour,  (W) 

fun 

unused  energy,  (Wh) 

pav 

char 

maximum  power  that  the  FC  together  with  the  battery 

H2avg 

average  hydrogen  level  of  the  tank,  (%) 

can  absorb,  (W) 

H2mjn 

minimum  hydrogen  level  of  the  tank,  (%) 

pav 

1  dis 

maximum  power  that  the  FC  together  with  the  battery 

Hfc 

electrolyzer  current  number  of  working  hours,  (h) 

can  deliver,  (W) 

Hfc 

FC  current  number  of  working  hours,  (h) 

pav 

1  H2,fc 

maximum  power  that  the  FC  can  deliver  in  the  next 

H 

:c 

working  hours  of  the  FC,  (h) 

hour,  (W) 

H 

ife 

fc 

FC  expected  life,  (h) 

pav 

1  H2,lz 

maximum  power  that  the  electrolyzer  can  store  in  the 

H 

ife 

fc 

FC  lifetime,  (h) 

next  hour,  (W) 

H 

ife 

Z 

electrolyzer  expected  life,  (h) 

Pb  at 

battery  output  power,  (W) 

H 

ife 

z 

electrolyzer  lifetime,  (h) 

pmax 

1  bat 

battery  maximum  power,  (W) 

U 

z 

working  hours  of  the  electrolyzer,  (h) 

P extra 

extra  power,  (W) 

HS 

hybrid  system 

Pfc 

FC  output  power,  (W) 

K 

warr 

FC  working  hours  defined  by  the  warranty,  (h) 

Pfcd 

FC  rated  power,  (W) 

f^bat.char 

battery  system  charge  constant,  (-) 

Plz 

electrolyzer  output  power,  (W) 

f^bat.des 

battery  system  discharge  constant,  (-) 

prtd 

electrolyzer  rated  power,  (W) 

f^H2,char 

hydrogen  system  charge  constant,  (-) 

Pnet 

net  power,  (W) 

f^H2,des 

hydrogen  system  discharge  constant,  (-) 

PV 

photovoltaic  solar  panel 

£h2 

level  of  the  hydrogen  tank,  (%) 

s 

state,  (-) 

Nbat 

number  of  batteries  used  throughout  the  lifetime  of 

soc 

state  of  charge 

^cycles 

the  HS,  (-) 

number  of  complete  cycles  of  the  battery,  (-) 

socavg 

SOCmin 

average  battery  SOC,  (%) 
minimum  SOC  allowed 

Nfc 

number  of  FC  used  throughout  the  lifetime  of  the  HS, 

t 

duration  of  the  battery  power  exchange,  (h) 

^cycles 

(-) 

maximum  charge/discharge  cycles  that  the  battery  can 
complete  during  its  life,  (-) 

At 

interval  time  between  two  simulation  steps,  (h) 

about  the  performance  of  off-grid  system  were  detailed  in  Refs. 

[12-14]. 

Furthermore,  due  to  the  imbalance  between  the  demanded 
power  and  the  power  generated  by  the  renewable  energy  system, 
an  ESC  is  necessary  for  the  optimal  energy  flow  of  the  energy 
sources  within  the  HS,  whose  main  goal  must  be  supply  the  load. 
Normally,  the  common  ESC  are  based  on  flowchart,  because  they 
are  easy  to  implement  and  really  effective.  Examples  of  ESC  based 
on  flowcharts  applied  to  HS  can  be  found  in  Refs.  [15-20  .  On  the 
other  hand,  ESC  not  based  on  flowcharts  have  been  also  used  [3,21 
23].  An  ESC  based  on  PI  controllers  was  used  in  Ref.  [3  .  A  simple 
ESC  based  on  the  available  energy  and  the  energy  demanded  by  the 
load  was  applied  in  Ref.  [21  .  An  ESC  based  on  mixed-integer  linear 
programming  was  detailed  in  Ref.  [22  .  A  genetic  algorithm  was 
developed  in  Ref.  23]  for  sizing  and  optimizing  the  energy  flow  of  a 
daily  load.  Nevertheless,  there  are  few  works  in  which  the  perfor¬ 
mance  of  the  ESC  have  been  studied  over  long  periods.  The  ESC  is 
evaluated  along  one  year  in  Refs.  18]  and  [23],  along  one  week  in 
Ref.  [24],  and  along  24  h  in  Ref.  [25]. 

This  work  presents  two  novel  ESC  for  off-grid  HS.  The  ESC  are 
applied  to  a  HS  composed  of  PV  panels,  WT,  and  battery  and 
hydrogen  system  (FC,  electrolyzer  and  hydrogen),  They  are  tested 


throughout  the  expected  life  of  the  HS  (25  years)  and  compared 
with  the  ESC  developed  by  the  authors  in  Ref.  [20]  with  the 
objective  of  reducing  the  utilization  costs  of  the  battery  and 
hydrogen  system.  The  first  ESC  tries  to  improve  the  generated  and 
stored  energy  of  the  HS,  whereas  the  second  one  is  focused  on 
minimizing  the  number  of  needed  elements  through  the  expected 
life  of  the  HS.  Thus,  once  concluded  the  introduction,  this  paper  is 
structured  as  follows.  The  configuration  of  the  off-grid  HS  used  for 
the  development  of  the  two  ESC  is  presented  in  Section  2.  Section  3 
details  the  two  novel  ESC,  and  develops  briefly  the  ESC  presented  in 
Ref.  20].  Section  4  presents  the  real  components  of  the  HS  used  to 
test  the  performance  of  the  ESC.  A  comparative  study  among  the 
three  ESC  is  carried  out  in  Section  5,  and  finally,  the  paper  ends  with 
the  conclusions  in  Section  6. 

2.  Off-grid  HS  configuration 

The  two  novel  ESC  can  be  applied  to  several  configurations  of  an 
off-grid  HS.  One  suitable  configuration  is  showed  in  Fig.  1,  which  is 
used  for  studying  and  comparing  the  ESC  presented  in  this  paper. 
According  to  this  configuration,  the  renewable  energy  system, 
composed  of  WT  and  PV  panels,  acts  as  primary  energy  source.  In 
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Turbine 


DC-bus 


Rectifier  + 
DC/DC 
Converter  + 
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^  Algoritlmi  y 


PV  System 


Load 


Hydrogen  tank 


Electrolyzer 


Fig.  1.  Off-grid  HS  configuration. 


both  cases,  the  WT  and  PV  panels  are  connected  to  the  DC  bus 
voltage  through  DC/DC  converters,  which  use  MPPT  controllers  to 
extract  the  maximum  available  power  from  the  wind  and  sun 
irradiance,  respectively. 

Otherwise,  the  battery  and  the  hydrogen  system  are  used  to 
support  the  renewable  energy  system  and  to  store  the  power 
generated  by  the  renewable  system  and  not  used  by  the  load.  Thus, 
if  the  demanded  load  is  higher  than  the  power  generated  by  the  WT 
and  PV  panels,  the  battery  or  the  FC  must  generate  this  power 
deficit.  On  the  contrary,  if  the  demanded  load  is  smaller  than  the 
renewable  power,  this  extra  power  is  used  by  the  electrolyzer  to 
generate  hydrogen  or  it  is  stored  in  the  battery.  In  all  the  cases,  in 
order  to  achieve  an  appropriate  energy  flow,  each  energy  source  is 
connected  to  the  DC  bus  through  a  controllable  DC/DC  power 
converter. 

3.  Energy  supervisory  controls  (ESC) 

As  it  was  commented  above,  the  power  generated  by  the 
renewable  energy  sources  does  not  correspond  sometimes  with  the 


demanded  load,  and  therefore  the  battery  or  the  hydrogen  system 
must  work.  The  difference  between  the  power  demanded  by  the 
load  and  the  power  generated  by  the  renewable  energy  sources  is 
defined  as  the  net  power,  which  can  be  positive  or  negative.  The 
ESC  has  to  decide  if  the  net  power  is  generated  or  absorbed  by  the 
battery  or  by  the  hydrogen  system  or  by  both,  and,  moreover,  to 
control  the  battery  SOC  and  the  hydrogen  level  between  optimum 
operating  margins. 

Fig.  2  presents  a  scheme  of  the  overall  configuration  used  for 
the  two  new  ESC.  The  main  difference  between  both  controls  is 
the  way  of  calculating  the  remaining  energy  of  the  battery 
(battery  charge  energy,  Ebat.char,  and  battery  discharge  energy, 
Tbat.dis)  and  hydrogen  system  (FC  energy,  Efc,  and  electrolyzer 
energy,  E\z).  Section  3.1  shows  the  common  blocks  for  both  ESC, 
and  Sections  3.2  and  3.3  detail  the  way  of  obtaining  the 
remaining  energies  in  the  two  ESC.  In  general,  both  ESC  try  to 
generate  or  absorb  the  net  power  with  the  battery  and  the 
hydrogen  system  by  applying  a  proportional  division  of  the  net 
power  between  them  taking  into  account  the  values  of  remaining 
energy,  as  shown  further  on. 


Fig.  2.  Scheme  of  the  overall  configuration  used  for  the  two  new  ESC. 
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Finally,  Section  3.4  describes  briefly  the  third  ESC,  already  pre¬ 
sented  in  Ref.  [20  .  This  third  ESC  based  on  flowchart  and  fuzzy 
logic  presents  a  different  approach.  It  was  designed  to  decrease  the 
utilization  costs  of  the  battery  and  hydrogen  system.  In  fact,  this 
ESC  and  the  other  two  novel  ESC,  in  which  the  utilization  costs  have 
not  been  taken  into  account,  are  compared  in  this  work  in  order  to 
show  the  differences  among  them. 


power  that  the  FC  together  with  the  battery  can  deliver,  whereas 
^char  t^ie  maximum  power  that  the  battery  and  the  electrolyzer 
can  absorb. 


pav  _  pav  I  pav 
1  dis  1  bat, dis  '  1  H2,fc 


pav  _  pav  i  pav 

1  char  1  bat, char  '  1  H2,lz 


3.1.  Development  of  the  ESC:  common  blocks 

According  to  Fig.  2,  the  common  blocks  of  the  new  ESC  are  the 
following  ones:  proportionality  constants,  ideal  power,  maximum 
available  power  calculation  and  flowchart.  Next,  all  these  blocks  are 
commented. 

3.1.1.  Proportionality  constants 

As  commented  above,  the  ESC  are  based  on  a  proportional  di¬ 
vision  of  the  net  power.  Thus,  in  both  controls  during  a  discharge, 
the  proportional  constants  are  obtained  as  follows: 


where  Pdat  charg,  Pdat  dis  are  battery  available  power  during  a  charge 
or  discharge  in  the  next  hour;  Pff2  fc  is  the  maximum  power  that  the 
FC  can  deliver  in  the  next  hour  (available  power  in  the  hydrogen 
tank);  and  Pff2  lz  is  the  maximum  power  that  the  electrolyzer  can 
store  in  the  tank  in  the  next  simulation  step.  Pdat  char,  Pdat  dis  can  be 
obtained  according  to  Eqs.  (10)  and  (11),  respectively  [26].  Pj^tchar 
is  the  minimum  value  between  the  maximum  power  and  the 
maximum  charge  power  of  the  battery,  which  depends  on  the  SOC. 
Pdat  dis  is  the  minimum  value  between  the  maximum  power  and  the 
maximum  power  that  the  battery  could  deliver  for  avoiding  a  SOC 
below  the  minimum  value  allowed  to  it  (SOCmin). 


/Ci 


bat, dis  — 


'bat, dis 


'bat, dis  “F  Efc 


(i) 


^bat,char  Him 


(cnom 
pnnax  ^bat 
bat  >  At 


x 


100 -SOC 
100 


(10) 


^H2,dis 


'bat, dis  "T  Efc 


(2) 


^bat,dis  mm 


(cnom 
pm  ax  ^bat 
bat  >  At 


SOC  -  SOC 


X 


min 


100 


(11) 


where  Kbat.dis  and  /<u2,dis  are  the  battery  and  the  hydrogen  system 
discharge  constants. 

In  the  same  way,  during  charging  process  the  proportional 
constants  are  the  following. 


^bat,char 

^bat,char 

(3) 

^bat,char  Fjz 

^H2,char 

 Elz 

(4) 

£bat,char  T-  E\z 

where  Kbarchar  and  I<m, char  are  the  battery  and  the  hydrogen  system 
discharge  constants. 


3.1.2.  Ideal  powers 

Once  obtained  the  proportional  constants,  the  ideal  battery 
power  and  the  ideal  hydrogen  system  power  (FC  and  electrolyzer 
power)  are  calculated.  These  powers  are  proportional  to  the  net 
power,  as  expressed  in  Eqs  (5)-(7),  and  they  would  be  the  power  to 
be  generated  by/stored  in  the  battery  and/or  hydrogen  system  if  the 
battery  SOC  and  the  hydrogen  tank  level  were  not  controlled. 


f  l^bat,dis  x  Aiet  with  Pne t  >  0 

\  l^bat,char  x  Enet  with  Pnet  <  0 


Pfc  —  ^H2,dis  x  Pnet  (6) 

P\z  =  ^H2,char  x  Enet  (7) 


3.1.3.  Maximum  available  powers  calculation 

The  maximum  available  powers  are  the  maximum  powers  of  the 
battery,  FC  and  electrolyzer,  which  are  defined  to  avoid  an  over 
charge  or  over  discharge  of  the  battery  and  the  hydrogen  tank 
during  the  next  simulation  step  (one  hour).  Furthermore,  the 
maximum  available  powers  in  a  discharge  and  in  a  charge  are 
defined  as  follows.  In  the  next  simulation  hour,  PjjV  is  the  maximum 


where  F£°m  represents  the  battery  nominal  capacity  in  energy  terms; 
P^x  is  the  battery  maximum  power  defined  by  the  manufacturer 
[27];  and  At  is  the  interval  time  between  two  simulation  steps. 

Regarding  the  hydrogen  system,  the  maximum  power  that  the 
FC  and  the  electrolyzer  can  deliver  or  absorb  during  the  next 
simulation  step  is  related  with  the  current  level  of  the  hydrogen 
tank.  Thereby,  the  FC  maximum  available  power  is  calculated  from 
the  minimum  value  between  the  FC  rated  power  (P£d)  and  the 
maximum  power  that  the  FC  could  generate  according  to  current 
level  of  the  hydrogen  tank  {Lm)- 

PH2,fc  =  min  (Pfcd  -  £1ow,H2  •  Vtherm '  Uf  ■  ??stacl< '  X 

(12) 


where  Fi0W,H2  is  the  hydrogen  lower  heating  value;  CAPh2  is  the 
tank  nominal  capacity  in  kilograms;  Vtherm,  Elf  and  ?7stack  are  specific 
parameters  of  the  FC,  which  depend  on  the  selected  FC;  Vtherm  is  the 
thermodynamic  efficiency;  Uf  is  the  FC  utilization  factor  [28];  and 
?7stack  is  the  FC  stack  efficiency,  which  is  calculated  as  the  ratio  be¬ 
tween  the  output  stack  voltage  and  the  standard-state  reversible 
voltage. 

With  respect  to  the  electrolyzer,  the  maximum  available  power 
that  can  be  absorbed  during  a  certain  hour  is  calculated  from  the 
minimum  value  between  the  electrolyzer  rated  power  (Pgd)  and  the 
power  that  could  produce  a  complete  filling  of  the  hydrogen  tank. 


m 2,iz  =  niin  (p[ztd .  B  ■  +  A 


CAPH2  /100-L 


H2 


At 


100 


(13) 


where  A  and  B  are  coefficients  of  the  selected  electrolyzer  model 
[29];  and  is  the  electrolyzer  rated  hydrogen. 


3.1.4.  Flowchart 

Finally,  because  of  the  battery  SOC  and  hydrogen  tank  level  must 
be  controlled,  a  controller  based  on  flowchart,  common  for  both 
ESC,  is  developed.  The  inputs  of  this  flowchart  are  the  ideal  power 
of  the  battery  and  hydrogen  system  (P£at,  P^andPj^),  the  net 
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power  (Pnet),  and  the  maximum  available  powers  during  a 
discharge  and  charge  (P^sand  Pa^ar).  outputs  are  the  battery,  FC 
and  electrolyzer  powers  (Pbat,  Pfc  and  P\z). 

Note  that,  by  using  the  maximum  available  powers,  the  battery 
SOC  and  the  hydrogen  level  can  be  controlled  and  the  power 
generated/absorbed  by  the  battery  and  by  the  hydrogen  system  can 
be  optimized.  Fig.  3  shows  the  flowchart  and  Table  1  presents  a 
summary  of  the  different  control  states  used  in  both  ESC.  During  a 
discharge  process  (Pne t  >  0),  the  flowchart  works  in  the  following 
way.  If  PgY  is  smaller  than  Pne t,  then  the  battery  and  FC  must 
generate  theirs  maximum  available  powers  (Pdatdis  andPf^fc)-  In 
this  state  (S  =  1),  the  HS  system  cannot  generate  the  power 
demanded  by  the  load,  and  therefore,  this  state  must  be  avoided. 
On  the  contrary,  if  P*[Y  is  higher  than  Pne t,  the  power  demanded  by 
load  can  be  generated.  Now,  if  the  battery  maximum  available 
power  is  higher  than  the  ideal  battery  power  (Pdatdis  >  Pdat)  but 
the  FC  maximum  power  is  smaller  than  the  FC  ideal  power  (P™2  fo  < 
PfC),  the  FC  must  deliver  its  maximum  power  and  the  battery  must 
generate  the  ideal  battery  power  plus  an  extra  power  (Pextra)  in 
order  to  provide  the  net  power.  This  is  named  as  state  2  (5  =  2). 
Otherwise,  if  the  previous  conditions  are  not  achieved,  the  opposite 
case  happens  (Pdatdis  <  P£at  and  Pg£  fc  >  Pj?c).  In  this  state  (S  =  3), 
the  battery  must  generate  its  maximum  available  power  and  the  FC 
has  to  generate  the  FC  ideal  power  plus  the  extra  power.  If  this 
latter  condition  neither  occurs,  then  the  maximum  available 
powers  are  higher  than  the  ideal  powers  (Pg^t  dis  >  P£at  and  P^2  fc  > 
PfC).  This  last  case  (S  =  4)  must  be  the  more  usual  state,  in  which  the 
battery  and  FC  way  deliver  their  ideal  powers. 

Finally,  if  Pnet  is  smaller  than  0,  the  flowchart  works  in  a  similar 
way  as  explained  above.  Thus,  the  battery  and  the  electrolyzer 
absorb  the  ideal  powers  only  if  the  following  two  conditions  are 


achieved:  < 


P*  I 
1  bat  I 


andP-<|Pi*z|,  (S  =  8). 


3.2.  ESC  \ 


Table  1 

Summary  of  the  control  states  used  in  the  flowchart  of  the  ESC  1  and  ESC  2. 
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calculates  the  remaining  energy  of  the  battery  (Ebat.char  and 
Pbat.dis)  and  hydrogen  system  (EfZ  and  Eiz)  from  the  battery  SOC 
and  the  hydrogen  tank.  In  this  case,  during  a  discharge,  the 
battery  available  energy  is  calculated  taking  into  account 
the  corresponding  energy  to  the  difference  between  the  current 
SOC  and  the  minimum  SOC  allowed  to  the  battery  (SOCmin).  In 
the  same  way,  during  a  charge,  the  energy  considered  by  the 
control  is  the  maximum  energy  that  the  battery  can  store  from 
the  current  SOC  to  a  complete  charge  of  the  battery.  For  this 
control,  all  the  energy  terms  are  denoted  as  ESC,1  in  the 
superscripts. 


fecs,i 

cbat,dis 

cECS,l 

cbat,char 


DODnom  p  /  c  r\  r'  q  \ 

Iqq  x  ^bat,nom  x  (SOC  —  SOCmjn) 

DODnom  r\r\  c  c'  \ 

TOO  X  ^at’nom  X  —  SOC) 


(14) 

(15) 


This  control  tries  to  improve  the  power  flow  of  the  HS 
avoiding  the  state  1  of  the  flowchart  and  to  store  most  of  the 
renewable  energy  in  the  hydrogen  tank  and/or  in  the  battery.  It 


where  Ebat.nom  is  the  battery  nominal  capacity  in  energy  terms;  and 
DODnom  represents  the  nominal  Depth  of  Discharge  (DOD)  rec¬ 
ommended  for  the  battery. 


Fig.  3.  Flowchart  control. 
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Regarding  the  hydrogen  system,  the  remaining  energies  are 
calculated  by  Eqs.  (16)  and  (17).  Eq.  (16)  represents  the  total  energy 
that  the  FC  could  generate  according  to  the  current  hydrogen  level. 
Eq.  (17)  defines  the  maximum  energy  that  the  electrolyzer  could 
store  in  the  hydrogen  tank  according  from  the  current  hydrogen 
level. 

Efc5’1  =  E1ow,H2  '  V  therm '  Ut  ■  V™?ck  ■  CAPH2  x  ^  (16) 

4CSI  =  B-CAPH2  +ACAPh2-(100100Lh2)  (17) 

where  V^ck  avera§e  stack  efficiency. 

Note  that,  with  this  ESC,  the  battery  SOC  and  the  hydrogen  level 
will  evolve  in  a  similar  way,  so  that  the  load  supply  could  be 
guaranteed  whenever  there  was  energy  available  stored  in  the 
battery  and  the  hydrogen  tank.  For  example,  during  a  positive  net 
power  (power  to  be  generated),  the  energy  source  with  more  en¬ 
ergy  stored  will  have  to  deliver  more  power  than  the  other  energy 
source  in  order  to  control  the  battery  SOC  and  the  hydrogen  tank 
level. 

3.3.  ESC  2 


divided  by  the  time  in  hours  (rtj^fe)  that  the  electrolyzer  takes  in 
completing  a  charge  cycle  of  the  hydrogen  tank. 
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3.4.  ESC  3 

The  third  ESC,  developed  in  detail  in  Ref.  [20],  attempts  to 
reduce  the  utilization  costs  of  the  battery  and  hydrogen  system.  It, 
denoted  as  ESC, 3  in  the  superscripts,  presents  a  control  scheme 
based  on  flowchart  and  fuzzy  logic.  The  flowchart  is  responsible  for 
deciding  the  power  to  be  generated  by/stored  in  the  battery  and/or 
hydrogen  system  while  trying  to  minimize  their  utilization  costs. 
The  reference  power  defined  by  the  flowchart  control  is  supervised 
by  a  fuzzy  logic  controller  in  order  to  keep  the  battery  SOC  and 
hydrogen  tank  level  between  certain  target  margins. 


This  control  tries  to  minimize  the  number  of  needed  elements 
(batteries,  FC  and  electrolyzers)  throughout  the  expected  life  of 
the  HS.  In  this  case,  the  remaining  energy  of  the  battery  and 
hydrogen  system  depend  on  the  total  quantity  of  energy  that  the 
battery  or  the  hydrogen  system  can  absorb  or  deliver  along  their 
whole  life,  which  decreases  as  the  elements  become  older. 
Thereby,  the  control  needs  to  know  parameters  such  as  the  cur¬ 
rent  number  of  complete  cycles  of  the  battery  (NCyCies)  and  the 
working  hours  of  the  FC  and  electrolyzer  (Hfc,  H\z).  Because  of  the 
battery  SOC  and  the  hydrogen  level  are  not  needed  to  calculate  the 
remaining  energy  terms,  the  flowchart  is  essential  to  control  these 
levels. 

The  battery  energy  terms,  denoted  for  this  control  with  ESC, 2  in 
the  superscripts,  are  obtained  from  Eqs.  (18)  and  (19).  In  this  case, 
the  considered  energy  corresponds  to  the  total  remaining  energy 
throughout  the  lifetime  of  the  battery.  Thus,  as  the  battery  is 
completing  cycles,  the  remaining  energy  is  decreasing  and  it  will  be 
used  less  than  at  the  beginning  of  its  life,  as  the  battery  gets  older. 
Moreover,  the  remaining  energy  is  considered  the  same  in  the 
charging  and  discharging  process. 


Sa'dis  =  05  ><  £bat, 
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where  N^les  are  the  maximum  cycles  of  charge  and  discharge  that 
the  battery  can  complete  during  its  life. 

The  FC  energy  is  obtained  from  the  product  of  the  total  energy 
that  the  FC  can  generate  from  the  hydrogen  stored  in  the  tank  -  the 
first  term  in  brackets  in  Eq.  (20)  -  and  the  difference  between  the 
FC  expected  life  in  working  hours  (H{jfe)  and  the  current  number  of 
working  hours  (Hfc),  which  is  divided  by  the  time  in  hours  that  the 
FC  takes  in  consuming  the  entire  hydrogen  tank  when  the  FC  is 
generating  its  rated  power. 

Similarly,  the  electrolyzer  energy  is  calculated  from  the  product 
of  the  maximum  quantity  of  energy  that  can  be  stored  in  the 
hydrogen  tank  -  the  first  term  in  brackets  in  Eq.  (21)  -  and  the 
difference  between  the  electrolyzer  expected  life  in  working  hours 
(H}*e)  and  the  current  number  of  working  hours  (Hiz),  which  is 


4.  Case  study:  off-grid  HS  located  in  Algeciras  (Spain) 

In  order  to  check  the  performance  of  both  ESC,  an  off-grid  sys¬ 
tem  located  in  Algeciras  (Cadiz)  Spain  have  been  considered.  The 
power  demanded  by  this  system  during  a  day  and  throughout  a 
year  is  shown  in  Fig.  4. 

The  consumption  curve  corresponds  to  a  load  in  which  the 
consumption  peak  occurs  from  12.00  to  16.00  pm  and  it  increases 
in  spring  and  summer  achieving  a  maximum  demanded  power 
around  1000  W. 

The  main  characteristics  of  the  battery,  hydrogen  system,  WT 
and  PV  system  are  developed  below.  The  sizing  of  these  compo¬ 
nents  has  been  performed  by  using  Simulink  Design  Optimization 
(SDO)  of  MATLAB  [30  ,  mainly  taking  into  account  that  the  gener¬ 
ated  power  must  satisfy  the  demanded  curve  and  extend  the  bat¬ 
tery  life. 
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Fig.  4.  Load  power  consumption  during:  a)  two  days;  and  b)  one  year. 
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4.2.  Battery 

The  selected  battery  for  the  off-grid  HS  is  a  660  Ah  lead-acid 
battery.  Its  commercial  name  is  BAE  SECURA  PVS,  model  PVS-660 
[27],  which  is  composed  by  6  cells  with  a  unit  voltage  of  2.23  V, 
so  that  the  battery  nominal  energy  (F[J°tm)  is  9500  Wh.  Otherwise, 
as  indicated  in  the  battery  datasheet,  its  maximum  charge  current 
is  330  A  with  a  nominal  charge  power  of  about  4.3  kW,  whose  value 
is  considered  during  the  discharge  process  (P{^x).  Regarding  the 
battery  life,  the  manufacturer  ensures  that  the  battery  can  complete 
3150  cycles  according  to  the  IEC  61427  cycles,  which  has  been 
considered  as  lifespan  criteria  (Afees). 

With  respect  to  the  battery  SOC,  the  manufacturer  recom¬ 
mends  avoid  deep  discharges  of  more  than  80%  (DODn0m)  repre¬ 
senting  a  battery  minimum  SOC  (SOCmin)  of  20%.  The  battery  SOC 
is  calculated  by  Eq.  (22).  Hourly,  it  is  obtained  by  means  of  an 
energy  balance  between  the  energy  provided  by  and  stored  in  the 
battery. 

SOC(%)  =  100- (l  -gnkE^at-f)  (22) 

where  Pbat  is  the  power  exchanged  by  the  battery;  and  t  is  the 
duration  of  the  power  exchange  (one  hour). 


4.2.  Hydrogen  system  (fuel  cell,  electrolyzer  and  hydrogen  tank) 


The  hydrogen  system  is  composed  of  the  following  commer¬ 
cially  available  components:  1)  The  FC  is  a  1.2  kW  polymer  elec¬ 
trolyte  membrane  (PEM)  FC  from  Heliocentris  [31],  model  NEXA 
1200;  2)  the  electrolyzer  is  the  model  HG-60  29]  from  Heliocentris, 
which  presents  a  nominal  power  consumption  of  480  W  with  a 
nominal  hydrogen  generation  of  60  standard  liters  h-1;  and  3)  the 
metal  hydride  hydrogen  tank  has  a  total  capacity  of  760  L. 

A  mass  balance  between  the  hydrogen  flow  generated  by  the 
electrolyzer  (qf^)  anc*  th e  hydrogen  flow  consumed  by  the  FC 
(qfg1)  is  carried  out  every  hour  in  order  to  obtain  the  hydrogen  tank 
level.  Compared  with  the  battery  SOC,  the  hydrogen  tank  can  be 
discharged  completely. 


LH2(%)  =  100 


(23) 


With  respect  to  the  expected  life  of  the  elements,  the  manu¬ 
facturer  does  not  provide  any  information  about  the  maximum 
working  hours.  The  manufacturer  gives  information  about  the 
degradation  rate.  In  the  case  of  the  FC,  the  output  voltage  is 
degraded  0.54  mV  per  each  working  hour,  with  an  initial  warranty 
of  1500  h.  Thus,  if  a  maximum  degradation  is  fixed  for  the  FC 
voltage,  the  FC  lifetime  can  be  calculated  by  Eq.  (24).  In  this  case, 
the  allowed  voltage  drop  has  been  fixed  to  a  value  which  produces  a 
very  low  output  voltage  but  enough  to  make  the  FC  DC/DC  con¬ 
verter  work  properly. 
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where  Hwarr  is  the  FC  working  hours  defined  by  the  warranty  [h]; 
V^ss  is  the  maximum  voltage  drop  allowed  in  the  FC  output  voltage 
[V];  and  Rfc  is  the  degradation  rate  of  the  FC  voltage  [V  hr1]. 

Similarly,  the  working  hours  deteriorate  the  performance  of  the 
electrolyzer.  Thus,  the  electrolyzer  will  need  to  absorb  more  energy 
at  the  end  of  its  life  than  at  the  beginning  in  order  to  produce  the 
same  hydrogen  flow.  Fixing  a  minimum  efficiency  in  the  electro¬ 
lyzer,  its  lifespan  in  hours  (H{j.fe)  can  be  calculated  as  follows. 


1 


(25) 


where  R\z  is  the  degradation  rate  of  the  electrolyzer  performance 
during  an  hour;  and  R^z  is  the  relationship  between  the  electrolyzer 
nominal  efficiency  at  beginning  of  its  life  and  the  minimum  effi¬ 
ciency  allowed. 


4.3.  Renewable  energy  system  (WT  and  PV  panels) 

In  the  proposed  HS,  the  renewable  energy  system  is  composed 
of  WT  and  PV  panels.  These  energy  sources  are  independent  of  the 
ESC  because  they  incorporate  DC/DC  converters  with  MPPT 
controller.  Both  energy  sources  generate  the  maximum  energy 
available  from  the  wind  and  sun  irradiance.  The  ESC  must  have 
information  every  hour  about  the  energy  generated  by  them. 

The  selected  WT  is  from  Bornay  [32]  with  a  nominal  power  of 
1.5  kW,  and  the  PV  system  consists  of  nine  EOPLLY  33]  panels 
connected  in  parallel,  with  a  total  peak  power  of  1.62  kW  p. 


5.  Simulation  results  and  discussion 

A  comparative  study  between  the  two  novel  ESC  presented  in 
this  work  (ESC  1  and  ESC  2)  and  the  third  ECS  based  on  reducing  the 
utilization  costs  (ESC  3)  is  carried  out  in  this  section.  The  controls 
are  tested  by  simulation  throughout  the  expected  life  (25  years)  of 
the  HS  showed  in  Fig.  1  with  the  components  described  in  Section 
4.  MATLAB-Simulink  was  used  in  the  modeling  and  performed 
simulations,  considering  a  sampling  time  of  one  hour. 

The  sun  irradiance  and  the  wind  speed  were  taken  throughout  a 
year  from  a  weather  station  located  in  Algeciras  (Cadiz),  Spain. 
These  data  (shown  in  Fig.  5)  and  the  demanded  load  showed  in 
Fig.  4  were  extrapolated  every  year  during  25  years. 

Fig.  6  (Fig.  7)  shows  the  power  exchanged  by  the  battery 
(hydrogen  system)  with  the  three  ESC.  Fig.  6a-c  (Fig.  7a-c) 
represent  the  battery  (hydrogen  system)  power  throughout  25 
years  in  the  three  ESC,  and  Fig.  6d  (Fig.  7d)  shows  a  detail  of  this 
power  during  two  weeks.  Fig.  8  shows  the  proportional  constants  of 
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Fig.  5.  a)  Wind  speed  (ms  1);  and  b)  sun  irradiance  (kW  m  2). 
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the  ESC  1  and  ESC  2.  In  both  controls,  the  battery  and  the  hydrogen 
power  are  proportional  to  their  constants  and,  in  the  case  of  ESC  1, 
they  are  directly  related  to  the  battery  SOC  and  to  the  hydrogen 
level  (Fig.  9). 

In  the  ESC  1,  the  battery  SOC  and  the  hydrogen  level  vary  in  a 
similar  way.  Thus,  during  the  discharges,  the  FC  constant  is  around 
0.68  and  the  battery  constant  around  0.32.  On  the  other  hand, 
during  the  charges,  it  occurs  the  following:  1)  if  the  battery  is 
charged  and  the  hydrogen  tank  is  not  filled  completely,  because  of 
the  battery  cannot  absorb  more  energy,  then  Kbat.char  =  1  and 
f<U2,char  =  0;  2)  if  the  hydrogen  tank  is  completely  charged  and  the 
battery  can  absorb  more  energy,  then  Kbat,char  =  0  and  /CH 2, char  =  1 ; 
and  3)  if  the  battery  and  the  hydrogen  tank  are  charged,  then 
f^bat.char  =  f^H2,char  =  0,  so  that  the  available  energy  from  the 
renewable  energy  sources  cannot  be  stored. 

Regarding  the  ESC  2,  because  the  proportional  constants 
depend  on  the  lifetime  of  every  component,  its  performance  can 
be  better  observed  throughout  the  simulation  of  25  years.  The 
lifetime  evolution  of  the  battery,  FC  and  electrolyzer  is  represented 
for  the  ESC  in  Fig.  10.  This  ESC  avoids  the  use  of  the  components  at 
the  end  of  their  life.  Thus,  Fig.  6b  shows  that,  from  around  the  year 
8.5  to  the  year  10  and  from  the  year  17  to  the  year  20,  the  battery 
is  less  used  and  the  FC  must  generate  more  power  than  the  rest  of 
the  time.  The  reason  for  this  fact  is  that,  during  the  time  periods 
indicated,  the  remaining  percentage  of  lifetime  of  battery,  and 
thus  the  remaining  available  energy,  is  very  small  compared  with 
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Fig.  6.  Battery  power:  a)  ESC  1  (25  years);  b)  ESC  2  (25  years);  c)  ESC  3  (25  years);  and 
d)  detail  of  two  weeks. 
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Fig.  7.  Hydrogen  system  power:  a)  ESC  1  (25  years);  b)  ESC  2  (25  years);  c)  ESC  3  (25 
years);  and  d)  detail  of  two  weeks. 


the  remaining  percentage  of  lifetime  of  the  FC  and  electrolyzer 
(Fig.  10).  Furthermore,  during  those  years,  the  battery  suffers 
lower  discharges  and  the  hydrogen  tank  reaches  very  low  levels 

(Fig.  9). 

It  can  be  observed  that,  with  the  ESC  3,  the  battery  has  to  work 
delivering  or  absorbing  small  values  of  power  (Fig.  6c),  and 
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Fig.  8.  Proportional  constants:  a)  ESC  1  (two  weeks);  and  b)  ESC  2  (25  years). 
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Fig.  9.  a)  Battery  SOC  with  ESC  1  (25  years);  b)  hydrogen  level  with  ESC  1  (25  years);  c)  battery  SOC  with  ESC  2  (25  years);  d)  hydrogen  level  with  ESC  2  (25  years);  e)  battery  SOC 
with  ESC  3  (25  years);  f)  hydrogen  level  with  ESC  3  (25  years);  g)  detail  of  battery  SOC  (two  weeks);  and  h)  detail  of  hydrogen  level  (two  weeks). 


therefore,  the  power  required  from  the  hydrogen  system  is  higher 
than  in  the  other  ESC  (Fig.  7c).  This  is  because  when  the  power  to  be 
absorbed  or  generated  is  small,  the  battery  utilization  cost  is  lower 
than  the  hydrogen  system  utilization  cost.  On  the  contrary,  if  the 
power  to  be  absorbed  or  generated  is  high,  the  cost  of  using  the 
hydrogen  system  is  cheaper  than  the  battery.  This  causes  that  the 
hydrogen  level  presents  deeper  discharges  than  the  battery  SOC,  as 
seen  in  Fig.  9e-f 

Finally,  the  following  parameters  were  calculated  and  studied  in 
order  to  compare  both  ESC: 

■  Number  of  batteries,  FC  and  electrolyzers  used  throughout  the 
lifetime  of  the  HS  (Nbat,  Nfc,  and  N\z)  and  total  number  of  com¬ 
ponents  used  (Nt). 

■  Unused  energy  (Fun),  which  is  defined  as  the  energy  available 
from  the  renewable  energy  sources  and  not  stored  in  the  battery 
and/or  in  the  hydrogen  tank. 

■  Total  time  percentage  that  the  ESC  have  been  working  in  the 
states  4  and  8  (P48),  time  percentage  in  the  state  1  (Pi)  and  time 
percentage  in  the  state  5  (P5). 

■  Average  and  minimum  values  of  the  battery  SOC  and  hydrogen 
level  (SOCmjn,  SOCaVg»  H2min,  and  H2ayg). 


■  Battery  efficiency  (^bat),  hydrogen  system  efficiency  (77H2)  and 
HS  efficiency  (?7hs),  which  were  calculated  by  Eqs.  (26)-(28). 

■  Total  acquisition  cost  (Act).  This  term  is  defined  as  the  sum  of  the 
acquisition  costs  of  the  battery,  FC  and  electrolyzer  along  25 
years  and  updated  in  every  replacement  by  means  of  the  annual 
interest  rate. 

■  Average  utilization  cost  (C^ed).  The  utilization  cost  is  calculated  by 
Eq.  (29),  in  which  the  utilization  costs  of  battery  (Cbat),  FC  (Cfc)  and 
electrolyzer  (Qz)  in  every  working  hour  are  added.  The  utilization 
costs  of  the  devices  were  calculated  as  defined  in  Ref.  [20]. 


(26) 


VH2 


^7lz,con '  ^fc,con 


P,zdt 


(27) 
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Fig.  10.  Lifetime  evolution  (25  years):  a)  battery:  b)  FC;  and  c)  electrolyzer. 


^HS 


^load^f 


[  (Pwt  +  Pp v)dt  +  Plow,H2  [  ^H2 dt  + 

JO  JO 


^bat^f 


(28) 


Qised  —  Q>at  +  Qc  +  Qz  (29) 

Table  2  shows  the  results  of  the  comparison  parameters  ob¬ 
tained  by  each  ESC.  According  to  these  results,  it  can  be  observed 
that,  as  expected,  the  ESC  2  needs  less  components  during  the 
lifetime  of  the  HS.  Nevertheless,  with  the  ESC  2,  the  time  per¬ 
centage  in  the  state  5,  and  thus  the  unused  energy,  are  higher  than 
with  the  ESC  1.  Respect  to  the  ESC  1,  the  working  time  in  the  states 
in  which  the  control  generates  the  ideal  powers  (states  4  and  8; 
p48  =  64.4%)  is  higher  than  with  the  ESC  2.  This  is  due  to  the  fact 
that,  in  the  ESC  1,  the  control  of  the  battery  SOC  and  hydrogen  level 
is  incorporated  in  the  own  control,  and  therefore  the  flowchart  is 
less  important.  On  the  other  hand,  with  respect  to  the  efficiency 
parameters,  the  battery  SOC  and  the  hydrogen  level  present  quite 
similar  values.  Hence,  the  HS  efficiency  reaches  a  value  slightly 
greater  than  50%  (a  little  greater  with  the  ESC  2),  and  the  battery 
SOC  is  never  below  the  20%.  Moreover,  the  minimum  total  acqui¬ 
sition  cost  is  achieved  with  the  ESC  2.  This  is  mainly  due  to  that  less 
replacement  are  required  with  this  ESC.  However  this  fact  does  not 
imply  that  the  average  utilization  cost  achieved  with  this  ESC  is  the 
lowest  of  all.  Finally,  as  expected,  the  lowest  average  utilization  cost 
is  obtained  with  the  ESC  3,  however,  this  ESC  presents  the  highest 
total  acquisition  cost  and  unused  energy. 


6.  Conclusions 

This  paper  has  presented  and  compared  two  novel  ESC  appli¬ 
cable  to  off-grid  HS  based  on  renewable  energy,  hydrogen  and 


Table  2 

Results  of  the  comparison  parameters  obtained  by  each  ESC. 


Parameter 

Value 

ESC  1 

ESC  2 

ESC  3 

Nbat 

3 

3 

5 

Nfc 

4 

3 

3 

Niz 

6 

5 

4 

Nt 

13 

11 

12 

Eu„  (kw  h) 

3650 

4334 

5443 

p48  (%) 

64.4 

54.5 

— 

p5  (%) 

28.6 

32.1 

— 

Pi  (%) 

0 

0 

— 

Vbat  (%) 

82.3 

81.2 

82.4 

VH2  (%) 

32.4 

32.5 

32.3 

Vhs  (%) 

53.2 

54.5 

53.2 

SOCmin  (%) 

21.2 

20.5 

20.6 

SOCavg  (%) 

94.9 

95.8 

82.94 

H2min  (%) 

0 

0 

0 

H2avg  (%) 

86.6 

87.5 

87.7 

Act 

52,716 

47,231 

50,780 

/-av 

'“used 

0.3467 

0.3143 

0.1523 

battery.  Furthermore,  these  ESC  were  also  compared  with  a  third 
ESC  (ESC  3),  whose  control  objective  was  to  reduce  the  utilization 
costs  of  the  energy  storage  devices. 

The  first  ESC  (ESC  1)  improved  the  energy  generated  by  the  HS 
and  stored  in  the  energy  storage  devices,  and  thus,  reducing  the 
non-used  energy.  On  the  contrary,  the  second  one  (ESC  2),  which 
considered  the  lifetime  of  the  HS  components,  reduced  the  number 
of  needed  components  along  the  expected  life  of  the  HS  and 
improved  the  HS  efficiency.  Moreover,  both  ESC  decided  the  battery 
power  and  the  hydrogen  system  power  (power  generated  by  the  FC 
or  power  absorbed  by  the  electrolyzer)  taking  into  account  the  net 
power  (difference  between  the  power  generated  by  the  renewable 
system  and  the  demanded  load),  the  battery  SOC,  the  hydrogen 
tank  level,  and  the  remaining  energy  of  the  battery  and  hydrogen 
system.  Depending  on  the  ESC,  these  energies  were  the  remaining 
energies  in  a  charge  or  discharge  cycle  (ESC  1)  or  the  total 
remaining  energies  which  could  be  delivered  or  absorbed  by  the 
components  during  their  lifetime  (ESC  2).  These  energy  parameters 
were  not  considered  in  the  ESC  3,  and  however,  the  utilization  costs 
of  the  energy  storage  devices  were  only  used  in  the  ESC  3. 

The  controls  were  tested  along  the  expected  lifetime  of  the  HS 
with  a  demanded  load  that  varied  depending  on  the  season  of  every 
year.  The  results  demonstrated  the  feasibility  of  these  ESC  for  off- 
grid  applications.  When  the  controls  were  compared,  the  ESC  2 
showed  better  results,  achieving  higher  HS  efficiency  than  the 
other  ESC,  and  reducing  the  number  of  needed  components  along 
the  expected  life  of  the  HS.  Nevertheless,  the  ESC  1  and  ESC  2  did 
not  improve  the  average  utilization  cost  obtained  with  the  ESC  3, 
since  the  utilization  costs  were  not  taken  into  account  in  these  ESC. 
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